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Abstract 


A refined zigzag theory is presented for laminated-composite and 
sandwich plates that includes the kinematics of first-order shear- 
deformation theory as its baseline. The theory > is variationally consistent 
and is derived from the virtual work principle. Novel piecewise-linear 
zigzag functions that provide a more realistic representation of the 
deformation states of transverse-shear-flexible plates than other similar 
theories are used. The formulation does not enforce full continuity > of the 
transverse shear stresses across the plate’s thickness, yet is robust. 
Transverse-shear correction factors are not required to yield accurate 
results. The theory > is devoid of the shortcomings inherent in the previous 
zigzag theories including shear-force inconsistency and difficulties in 
simulating clamped boundary conditions, which have greatly limited the 
accuracy of these theories. This new theory requires only C° -continuous 
kinematic approximations and is perfectly suited for developing 
computationally efficient finite elements. The theory should be useful for- 
obtaining relatively efficient, accurate estimates of structural response 
needed to design high-performance load-bearing aerospace structures. 


Nomenclature 

a, b 
2 h 

2 h (k) 


N 

S m 

S u > s* 
c° 


lateral dimensions of a rectangular plate 
total plate (laminate) thickness 
thickness of the /cth material layer (lamina) 

reference plate-coordinate axes positioned in the middle plane of the 
laminate 

thickness coordinate axis 

number of material layers (laminae) through the laminate thickness 
reference middle plane of the laminate 

parts of the cylindrical edge surface of the laminate where displacements 
and tractions are prescribed, respectively 

denotes a continuous function whose first-order derivative is 
discontinuous 
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c c 

U ’ '~'(7 

s,n 

q 


intersections of the cylindrical edge surfaces ( S u , 5^ ) with the middle 
surface S m where displacements and traction resultants are prescribed, 
respectively 

unit outward tangential and normal vectors to the mid-plane boundary 
(see Figure 1) 

applied transverse pressure [force/length 2 ] (see Figure 1) 


(T 13 T 2 ,T M ) 

(u[ k) ,u^\u z ) 


prescribed in-plane and transverse shear tractions (see Eq. (14)) 

in-plane and transverse components of the displacement vector in the /cth 
material layer (see Eqs. (1)) 


(u,v, w, 0 V 0 2 , y/ x , y/ 2 ) 


kinematic variables of the refined zigzag plate theory (see Eqs. (1)) 


(N l , N 2 , N n ) 


membrane stress resultants (see Eqs. (16)) 


(M t , M 2 ,M n ) 
(M'.mJ.M'.M') 

( Q„Qi ) 

€ (a = 1 , 2 ) 

P? (a = 1,2) 


"(*) 


(“(*)’ W 


rV,rt\rV) 


bending and twisting stress resultants (see Eqs. (17)) 

bending and twisting stress resultants due to zigzag kinematics (see Eqs. 

(17)) 

transverse shear stress resultants (see Eqs. (18)) 

transverse shear stress resultants due to zigzag kinematics (see Eqs. (18)) 
zigzag functions (see Eqs. (1)) 

derivatives of zigzag functions with respect to the thickness coordinate 
(see Eq. (2.2)) 

dimensionless thickness coordinates of the /cth layer (lamina) (see Eq. 

(5)) 

thickness coordinate of the interface between the /cth and ( k + 1) th 
layers (see Figure 1) 

dimensionless in-plane displacements along the interface between the /cth 
and ( k + 1) th layers (see Figure 2) 

strains in the /cth layer (see Eqs. (3)) 
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r a ( x i’X 2 ) (a = 1,2) 


average shear strains (see Eq. (2.1) and Eq. (9)) 


n a (« = i,2) 

¥ai x v x i) (« = 1,2) 

icr (k) a ik) r (k) r (k) r (i)N > 

\^11 ’ u 22 ’M2 > l 2 z ’ Mz ) 

E fk) 


Gf 




transverse-shear strain measures (see Eqs. (11)) 
zigzag amplitude functions (see Eqs. (1)) 
stresses in the Mi layer (see Eqs. (3)) 

Young's moduli of the Mi layer (see Table 1) 
shear moduli of the Mi layer (see Table 1) 
Poisson ratios of the Mi layer (see Table 1) 


A ,, , B , D : , and G„ 


constitutive stiffness coefficients (see Eqs. (22)) 



in-plane and transverse shear elastic stiffness coefficients for the /cth 
layer (see Eqs. (3)) 


(a = 1,2) 

£ 

<>(*),« 

k 2 


weighted-average, laminate-dependent transverse shear constants (see 
Eqs. (12.3)) 

variational operator (see eq. (14)) 
partial differentiation 

shear correction factor for First-order Shear Deformation Theory (FSDT) 
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1. Introduction 


High-performance and lightweight characteristics of advanced composite materials have spurred a 
wider range of applicability of these materials in military and civilian aircraft, aerospace vehicles, and 
naval and civil structures. To realize the full potential of composite structures, further advances in 
structural design and analysis methods are necessary. In particular, development of cost-effective and 
reliable laminated-composite structures as the primary load-bearing components of a vehicle requires 
further advances in stress analysis and failure prediction methodologies. 

A wide variety of modem civilian and military aircraft employ relatively thick laminated composites, 
with one hundred or more layers, in the primary load-bearing structures. Such structures can exhibit 
pronounced transverse shear deformation and, under certain conditions, design-critical thickness-stretch 
deformations. Fail-safe design of these structures requires accurate stress-analysis methods, particularly 
for regions of stress concentration. Computationally efficient analytical models based on beam, plate and 
shell assumptions that account for transverse shear and thickness-stretch deformations have recently been 
addressed in [1-3]. To achieve accurate computational models, three-dimensional finite element analyses 
are often preferred over beam, plate, and shell models that are based on First-order Shear-Deformation 
Theories (FSDT). This preference is because the latter tend to underestimate the normal stresses, 
particularly in highly heterogeneous and thick composite and sandwich laminates [4-7], For composite 
laminates with hundreds of layers, however, three-dimensional modeling becomes prohibitively 
expensive, especially for nonlinear and progressive failure analyses. To realize improved response 
predictions based on beam, plate, and shell assumptions, a variety of refined theories have been 
developed, e.g., [8-11]. Many of these refined theories have significant flaws in their theoretical 
foundation and predictive capabilities and, for these reasons, have not found general acceptance in 
practical applications. 

One class of refined theories that has emerged as practical for engineering applications is known as 
zigzag theories [12-24]. This class of theories employs a zigzag-like distribution for the in-plane 
displacements through the laminate thickness, while ensuring a fixed number of kinematic variables 
regardless of the number of material layers. It has been shown in [12-20] that zigzag theories provide 
sufficiently accurate response predictions for relatively thick laminated-composite and sandwich 
structures, including those for normal strains and stresses. Furthermore, these theories often yield the 
response predictions comparable to those that can be obtained from other layer-wise and higher-order 
theories that are more computationally intensive. To make a zigzag theory practical for a large-scale 
analysis and engineering design, the analytic framework of the theory must be well suited for an efficient 
finite element approximation. 

Recently, Tessler et al. [25] elucidated several serious flaws in the most notable zigzag theories (see 
the discussion in Section 3), and proposed a refined zigzag methodology that eliminates these flaws in an 
original and theoretically consistent manner. This new methodology is general enough to provide a means 
for advancing this class of theories and is expected to be significantly more attractive for engineering 
practice, including applications of these theories in large-scale finite element analyses. 

The aim of this report is to present a new zigzag theory for laminated-composite and sandwich plate 
structures that may exhibit a high degree of transverse shear flexibility, anisotropy, and heterogeneity. 
The range of applicability of the theory spans thin to relatively thick plates. This new theory is built upon 
the basic ideas of a refined zigzag theory for beams presented in [25]. Herein, FSDT is augmented with an 
improved zigzag kinematic field that involves a novel ^-continuous (across lamina interfaces) 
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representation of the in-plane displacements. The kinematic field is independent of the number of material 
layers and does not require enforcement of transverse-shear-stress continuity to yield accurate results. 
Unlike other similar theories (e.g., [12,21]), the zigzag contribution to the in-plane displacement field is 
physically realistic, is zero-valued at the top and bottom plate surfaces, and accounts for the shear 
deformation of every lamina in a consistent way. As a result, transverse-shear correction factors are not 
needed. Additionally, the plate equilibrium equations, constitutive equations, boundary conditions, and 
strain-displacement relations are consistently derived from the virtual work principle. 

One major benefit arising from the analytical form of this new theory is its ideal suitability to finite 
element modeling, where the kinematic-variable approximations need not exceed C° continuity. The 
implication of this attribute is that computationally efficient finite elements can be developed and used in 
general-purpose finite element codes. This benefit will enable an efficient use of accurate zigzag 
approximations in large-scale analyses to facilitate the development of robust designs of high- 
performance aerospace vehicles. 

In the remainder of the paper, the theoretical foundation of the new theory and its quantitative 
assessment are detailed. The zigzag kinematic assumptions, strain-displacement equations, and 
constitutive lamina relations are presented in Section 2. The major theoretical anomalies associated with 
the previous zigzag theories are briefly discussed in Section 3. A set of unique zigzag functions is then 
introduced and their mathematical structure is described in Section 4. The plate equilibrium equations and 
their associated boundary conditions, derived from the virtual work principle, are presented in Section 5. 
In Section 6, an extensive quantitative assessment of the theory is carried out, using closed-form solutions 
for simply supported and cantilevered plates made of laminated-composite and sandwich material 
systems. Some of the example problems represent significant challenges for any approximate theory. It is 
demonstrated that this new zigzag theory eliminates a major flaw of other similar theories; that is, the 
theory enables accurate modeling of the clamped boundary condition while adhering strictly to the 
variationally required boundary conditions. Finally, several concluding remarks emphasizing the merits of 
the new theory are presented in Section 7. 


2. Kinematics and formulation 


Consider a laminated plate of uniform thickness 2 h with N perfectly bonded orthotropic layers (or 
lamina) as shown in Figure 1. Points of the plate are located by the orthogonal Cartesian coordinates 
(x 1 ,x 2 ,z) . The ordered pair (x u x 2 ) e S m denotes the in-plane coordinates, where S m represents the set 
of points given by the intersection of the plate with the plane z = 0 , referred to herein as the middle 
reference plane (or midplane). The symbol z e[-h, h] is the through-the-thickness coordinate, with 
z = 0 identifying the midplane. The plate is subjected to a normal-pressure loading, q(x l ,x 2 ) , attributed 
to the midplane, S m , that is defined as positive in the positive z direction. In addition, a traction vector, 
( 7] , T 2 ,T_) , is prescribed on S a c S , where S denotes the total cylindrical-edge surface. On the 
remaining part of the edge surface, S u c= S , displacement restraints are imposed (or prescribed). The 
sections of the plate edge are related by S a u S u = S and S a C\ S u = 0 . Moreover, the curves 
C a = S a C\ S m and C u = S u n S m define the two parts of the total perimeter C = C a u C u surrounding 
the midplane region, S m . Finally, it is presumed that the lamina constitutive properties may differ 
appreciably from lamina to lamina, the plate deformations result in small strains, and that body and 
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inertial forces are negligible. 

The orthogonal components of the displacement vector, corresponding to material points of the plate 
(or laminate), are expressed as 


u[ k) (x 1? x 2 ,z) = u(x l ,x 2 ) + z 6* 1 (x 1 ,x 2 ) + ^ </l> (z) y/ l {x x ,x 1 ) 
U ( 2 \x { ,X 2 ,z) = v(x l5 x 2 ) + z 6* 2 (x 1 ,x 2 ) + (j) 2 k \z ) ^ 2 (Xj,X 2 ) 
u z (x l ,x 2 ,z) = w(x x ,x 2 ) 


where the in-plane displacement components u[ k 1 and u 2 k) are comprised of constant, linear, and zigzag 
variations through the thickness. The zigzag variations are C° -continuous functions with discontinuous 
thickness-direction derivatives along the lamina interfaces. The superscript ( k ) is used to indicate 
quantities corresponding to the Mi lamina, whereas the subscript (k) defines quantities corresponding to 
the interface between the k and (k+1) laminae. Thus, the kth lamina thickness is defined in the range 
z e [z (jtl) ,z (i) ] ( k = 1,...,^) (see Figure 1). The transverse displacement u z is assumed to be constant 

through the thickness and is independent of constitutive properties of the kth lamina; hence the superscript 
(k) does not appear in its definition. 



The kinematic variables in Eqs. (1) can be interpreted as follows. For homogeneous plates, the zigzag 
functions </>^ k) (a = 1,2) vanish identically and Eqs. (1) yield the kinematics of FSDT. For this 
degenerate case, u and v represent the midplane displacements along the coordinate directions x x and 
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x 2 , respectively; 0 l and 0 2 represent average bending rotations of the transverse normal about the 
positive x 2 and the negative x, directions, respectively; and w is the transverse deflection. For more 
precise definitions of the kinematic variables within FSDT refer, for example, to Ref. [1], The symbols 
(a = 1,2) denote through-the-thickness piecewise-linear zigzag functions associated with 
nonhomogeneous plates, yet to be defined. The y/ a = i// a (x l ,x 2 ) (cc = 1,2) functions represent the 
spatial amplitudes of the zigzag displacements and, together with the other five kinematic variables, are 
the unknowns in the analysis. The zigzag displacements 0P 1 y/ u (a = 1,2) may be regarded as 
corrections to the in-plane displacements associated with laminate heterogeneity. 

Consistent with the kinematic assumptions in Eqs. (1), the theory accounts for transverse shear 
deformation. (Transverse normal deformations are neglected in the kinematics; however, their inclusion 
may be possible following, for example, [1]). Correspondingly, the in-plane and transverse shear strains 
are 


s u ^ — u ,\ 'f z @\\ + 0\ ^V\,\ (2-1) 

S 22 ^ = V 2 + zQ 2 9 +^2 ^ 2 2 ( 2 . 2 ) 

Y\2 = W ,2 + V ,1 + z (^l,2 + ) "b + 02 ^2,1 (2-3) 

Yaz =Ya + Pa V« (a = 1,2) (2.4) 


d(-) 


where, henceforward, (•) a = denotes a partial derivative with respect to the midplane coordinate, 

’ dx a 

x a ( a = 1,2) . Also, the following notation is introduced 


Ya= w ,a+ 0 a (O' = 1,2) (2.5) 

(a = 1,2) ( 2 . 6 ) 

where the shear angles y a are uniform through the total laminate thickness, and 1 are piecewise constant 
functions that are uniform through the thickness of each individual lamina. 

The generalized Flooke's law for the Mi orthotropic lamina, whose principal material directions are arbitrary 
with respect to the midplane reference coordinates, (x, ,x 2 ) e S m , is written as 
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'°u 

(*) 

'C„ 

C 12 

Q« 

0 

0 ' 

(*) 


<J 22 


C 12 

c 22 

C 26 

0 

0 


8 22 

T \2 

= 

Q« 

c 26 

q 6 

0 

0 

« 

7X2 

*2z 


0 

0 

0 

Q 22 

Q\2 


72z 

. 7 lz . 


0 

0 

0 

Q\2 

Qn_ 


7\z . 


( 3 ) 


where Cf ] ( i,j = 1,2,6) and Q^j ( p,q = 1,2) are the transformed elastic stiffness coefficients referred to the 


(x 1 ,X 2 ,z) coordinate system and relative to the plane-stress condition that ignores the transverse-normal stress. The 

expressions for these coefficients in terms of the elastic moduli corresponding to the material coordinates can be 
found, e.g., in [10]. 


3. Anomalies of previous zigzag theories 


The ‘linear zigzag model’, that was developed for plate bending problems by Di Sciuva [12] and later 
examined further in [13-24], employs only five kinematic variables. In reference to Eqs. (1), these 
variables are u , v, w , y/\ an d y/ 2 . In contrast to Eqs. (1), Di Sciuva’ s zigzag-amplitude variables, (//, 

and y/ 2 , both contribute to the individual in-plane displacements, i.e., u\ k) = j\(u, vv , , , i// 2 ) and 

l{ i k) = / 2 (v, w , , \j /\ , ¥ 2 ) > thus resulting in a system of plate equilibrium equations that is coupled with 
respect to i//| and i// 2 (refer to Ref. [12] for the precise form of the displacement expansions of Di 

Sciuva’s theory). Di Sciuva determines the appropriate zigzag functions by enforcing the transverse 
shear stresses to be continuous along the adjacent lamina interfaces and, in addition, by requiring that the 
zigzag functions vanish in a single layer that is selected a priori, labeled herein as a ‘fixed layer.’ As a 
result of the above stated assumptions and transverse shear stress-constraint conditions, Di Sciuva’s 
theory attains the following characteristics: (a) the transverse shear stresses are uniform through the 
thickness and they correspond to those in the ‘fixed layer’; (b) the laminate transverse shear stiffness is 
governed by the transverse shear moduli of the ‘fixed layer’ alone; (c) the transverse shear strains and 
stresses erroneously vanish along fully clamped edges; and (d) the integral of the transverse shear stress 
across the laminate thickness does not correspond to the shear force resultant obtained from the 
corresponding plate equilibrium equations - a physical inconsistency (or anomaly) in the definition of the 
shear force. 

To remove the flaws associated with Di Sciuva’s ‘linear zigzag model’, Tessler et al. [25] introduced a 
refined zigzag theory for laminated-composite and sandwich beams in which: (1) a novel zigzag function 
is used to produce non-vanishing zigzag displacements in every lamina, thus removing the shear stiffness 
bias associated with the ‘fixed layer’ approach, and (2) the equilibrium of transverse shear stresses along 
adjacent lamina interfaces is fulfilled only in an average sense. The resulting theory is devoid of the 
aforementioned flaws of the previous zigzag theories and has been shown to demonstrate consistently 
superior results. In what follows, an extension of the methodology presented in [25] to plates is described 
in detail. 
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4. Refined zigzag functions and transverse shear constitutive relations 


The refined zigzag functions of the present theory are defined by using functions that are piecewise 
linear and C° continuous through the laminate thickness. For convenience, dimensionless lamina- 
interface displacements u (i) and v (l) (i = 0, 1,..., /V) are used to define fi '" 1 and (f ) 2 k 1 , respectively (see 
Figure 2 depicting the notation for a three-layered laminate). 

Thus, for the Mi lamina located in the range [z (Jt _ I)5 z (i) ] , (f>\ k ' and (/>-! 1 are given as 


^*) = l(l_|(*))„ (i _ |) + l(l + ^) %) (4) 

4 m =Ki-f m )v (i -, ) + Ki+# w >v w 

where 

=[(z- Z (k ^)/h (k) - 1] e [-1,1] (k = 1 (5) 


with the first lamina beginning at z {0) = — h, the last (ATh) lamina ending at z ( V) = h , and the 
A:th lamina ending at z {k} = z (k 1( + 2/; (/ 1 , where 2h {k) denotes the Arth lamina thickness. 



(a) Layer notation (b) Zigzag function <f)\ k 1 (c) Zigzag function <f^ k ] 

Figure 2. Notation for a three-layered laminate and (f)\ k * and (f)\ k) zigzag functions defined in terms 
of interfacial displacements, u (i) and v (i) (i = 0,1 . 
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Evaluating Eqs. (4) at the laminae interfaces gives rise to the definitions of the interfacial displacements 


= fC(f m = - 1), ««*) = = i) 

V (M> = — 1), v m = #’(£<*> = 1) 

where the interfacial displacements at the bottom 
identically, i.e. 

M (o, = M <v) = v m = V (N) = 0 ( 6 - 2 ) 

Substituting Eqs. (4) into Eqs. (2.6) results in the piecewise-constant functions {3 ( J C) given by 


(6.1) 
(k = 

and top plate surfaces are set herein to vanish 


[A 1 ' 1 ] _ 1 \ U U<) u (k- 1)1 . , 

{/?:>'} 2A«’{v m -v„_ l) J 

Because the zigzag functions are zero-valued on the top and bottom surfaces, as defined by Eqs. (6), 
through-the-thickness integrals of the slope functions (3 {k) ( a = 1,2) vanish identically, i.e. 


f N 



(8) 


U=i J 

Integrating Eqs. (2.4) across the laminate thickness and normalizing the result by the total laminate 
thickness, reveals that 



(9) 


Thus, y a (<2 = 1,2) represent the average transverse shear strains, coinciding with the common 
representation of the transverse shear strains used in FSDT. Also, Eqs. (9) indicate that the zigzag 
amplitude variables, y/ a {a = 1,2) , do not contribute to the average transverse shear strains. 

The u (k) and v ((t) interfacial displacements are obtained from Eqs. (7) in terms of 1 
{a = 1,2; k = 1,..., N) , i.e. 


■to I U w 




*(k- 1) 


v (k- 1) 


{k = \,...,N) 


( 10 . 1 ) 
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or, alternatively, 


u 


u-> 


V 


(*) 





( 10 . 2 ) 


with Eqs. (6.2) defining zero-valued top- and bottom-surface interfacial displacements. 

Following the approach in [25], the /3 { a k) functions are determined by first casting the transverse shear 
strains, Eqs. (2.4), in terms of the transverse-shear strain measures, rj a = y a —y/ a (a = 1,2) , and the 
zigzag amplitude functions, y/ a ( a = 1,2) , as 


Y\z 


(*) 


Yiz 


Eijjo + A 1 * 1 ) 

W 0 


0 

(i+A w ) 



(ii) 


The ij a strain measures are set to vanish explicitly in Di Sciuva’s theory [12] and enforced to vanish by 
way of penalty constraints in AverilTs theory [21], thus equating y a to y/ a . Presently, no such 
constraints are imposed on these strain measures. 

The transverse shear stresses using Eqs. (3) and Eqs. (1 1) are given as 


(k) r 


1 2 r _ 


011 012 
012 022 


l(h 


7i 

Jh , 


+ 


(i+A w ) o 
o d + /?f , )_ 


I ¥\ 

Wi, 


or, alternatively, they can be expressed as 


(k) r 


l lz 


1 2z , 


011 012 
012 022 


lW 


\vi 


( 12 . 1 ) 


+a ( J ) (i+A, 1 ") 


(*)>! 


I0i ( | ) /0i ( f ) 


'W\ + 


0«( 1 + A W ) 


Qvi / 0 (i) 


22 




( 12 . 2 ) 


In this form of the transverse-shear constitutive relations, the stress vector associated with the //^ strain 
measures is independent of the zigzag functions. The second and third stress vectors include, as their 
normalization factors, the coefficients 0^(1 + /?^) (a = 1,2) that are dependent on the zigzag 

functions through 1 . Elerein, these normalization factors are set to be constant quantities, denoted as 
G a ( a = 1,2) , thus imposing constraint conditions on the distribution of the zigzag functions. 
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These constraints give rise to the expressions for 


P\ 


(*) 

a 





G, 


0-22 


1 

1 


(12.3) 


The G x and G 1 constants are obtained by integrating Eqs. (12.3) through the laminate thickness while 
making use of Eqs. (8), resulting in 




h dz 
2 h^-h~Q^ 




h dz 
h~Q^) 


22 J 


Y 1 ' 


f 1 * h m ^ 

-1' 

J 


h^O (k) 


y 1 

> = < 

( 1 " h {k) ^ 

-1 

J j 


v n k= 1 1^22 J 



(12.4) 


where it is seen that G x and G 1 are weighted-average transverse-shear stiffness coefficients of their 
respective lamina-level coefficients, and Q 22 ■ 


Replacing Eqs. (12.3) into Eqs. (12. 1), results in the transverse-shear constitutive relations of the form 


yir 

“011 

012 

(*) 

< 

l r 2zj 

_012 

022 J 



f 


Y\+Y\ 

Y2+Y2 


Gl_i 

V0lY 7 
f r \ 
(*) 


A 


022' 


J 


(13.1) 


where, with the use of Eqs. (12.4), the dimensionless stiffness ratios are given as 


( nG) 


nW 

"daa V 


Odd [ h dz 




(a = 1 , 2 ) 


(13.2) 


These ratios are, in general, piecewise-constant through the laminate thickness; however, for 
homogeneous plates, they are unit-valued. Thus, for homogeneous plates, the zigzag transverse-shear 
contributions vanish, in which case Eqs. (13.1) become identical to the corresponding relations of FSDT. 


12 



The <f> ( a k) ( a = 1,2) zigzag functions are determined by substituting Eqs. (10.2) and Eqs. (12.3) into 
Eqs. (4), while making use of Eqs. (5), resulting in 


£ } =(z + h) 


g (1) 

\^aa 


-1 


(* = 1) 


=( z + h) I 


0, 


-1 


lx 


M) 


e; 


(M) 

aa 


Q w 


[k = 2,...,N) 


aa J 


(13.3) 


Z ^ l Z (k-l)’ Z (k)]’ Z ( 0 ) h’ Z (k) Z (k- 1 ) + 2/j (A' 1,...,^V , a 1,2) 

It is seen that the zigzag functions are independent of the state of deformation and are represented by C°- 
continuous, piecewise linear functions of the thickness coordinate. These functions possess slope 
discontinuity along lamina interfaces and vary linearly through the thickness of each lamina. Within each 
lamina, the transverse shear properties are represented by the stiffness ratios defined in Eqs. (13.2). For 
homogeneous plates, the zigzag functions vanish identically. 

The zigzag amplitudes, y/ a (a = 1,2), are vector functions of the actual response due to the applied 

loading, and they provide the proper scaling of the zigzag functions, thus controlling the total zigzag 
contribution to the in-plane displacements. The two zigzag amplitude functions and the remaining five 
kinematic variables constitute a set of seven kinematic variables associated with this refined zigzag plate 
theory. 


5. Equilibrium equations, boundary conditions, and constitutive relations 


The plate equilibrium equations and boundary conditions are derived from the virtual work principle 
which, neglecting body forces and assuming zero shear tractions on the top and bottom bounding plate 
surfaces, may be written as 


o = | } + 4* wf + ax' + * ds 

S m ~h 

h 

-\(dSw) d S-\\ \_T x Su[ k) + T 2 Su ( 2 k) + T z Su {k) ~\dsdz 

s w c, -h 

where 8 is the variational operator; all other symbols have been defined in Section 2. 


(14) 


f3 



Substituting Eqs. (1) and (2) into Eq. (14), and integrating across the plate thickness yields the two- 
dimensional statement of virtual work 

0 = f \N x 8u x + N 2 8v 2 + N X2 (8u 2 + 8v x ) 

+M X 80 XX + M 2 86 22 + M X2 (80 X2 + S0 21 ) 

+Q l (8w l + 89 x ) + Q 2 (8 w 2 + S0 2 ) - qSw]dS 

~\ c l N ln 8u + N ln 8v + QJw + M Xn 80 y + M. ln 80 x ]ds ( 1 5) 

+ [ \Mf8ys j a + M 2 8ys 2 2 + Mf 2 8i// X2 + M 2] 8i// 2 , + Q x 8y/ x + Q 2 8i// 2 ]dS 

-L [ M In SW\ + M 2n 8 iff 2 \ ds 

where the last two lines correspond to the zigzag kinematics contributions. In Eq. (15), the membrane 
stress resultants and conjugate strain measures are 


Nl - {ZV„ N 2 , N i2 } = J ' {<>, <rg\ rg>}& 

1. V ,2. ",2+ V ,l} 


Likewise, the bending stress resultants and conjugate strain measures are 


M T h = \M X , Mf, M 2 , Mt, M n , M* 2 , Mt_ x 


= |_ /! {zcr 1 ( f ) , (f{ k) off, za { 22 \ k) (J 22 , zr^ ,(j)[ k) t[ 1\ (fc k) r[^dz 

e b = {^1,1 5 ¥\,\ 5 @2,2 5 ¥2,2 5 @\,2 + ^2,1 2 ¥\,li Vix } 
and the transverse shear stress resultants and conjugate strain measures are 

Qj - { 02 , Qt a. Qt}= J>g>, rg\ r£\ r(f}& 

e i = { W .2 + ^2.K'2. M ',l+A.( / l} 

The force and moment resultants due to the prescribed tractions have the form 


N u ,N 2n ,Q zn ,M u ,M 2n ,Mf n ,Mi n 


(16.1) 


(16.2) 


(17.1) 


(17.2) 


(18.1) 


(18.2) 


= \ h h {T x ,T 2 ,T : ,zT x ,zT 2 , </\ k % 4 k) T 2 }dz 
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Integrating Eq. (15) by parts results in seven equilibrium equations and associated boundary conditions. 
The equilibrium equations are 


du : 

+ Sf X2 2 - 0 


Sv : 

N x2 j + N 2 2 = 0 


Sw : 

0i,i + 02,2 + q = o 


SO, : 

M U +M n, 2 ~Q\ =0 

(20) 

S0 2 : 

M X2 \ + ^2,2 ~ 02 = ^ 


Sy / X : 

K + M ki-Qt = o 


Sy / 2 : 

Mfn+A^f, -o:’=o 



The kinematic and force boundary conditions are given by 


u = u on C u or N x n x + N n n 2 = N Xn on C ff (21.1) 

v = v on C u or N n n, + N 2 n 2 = N 2n on C a (21.2) 

w = w on C u or Q 1 n x + Q 2 n 2 = Q zn on C CT (21.3) 

0 X = 0 X on C u or M x n x + M X1 n 2 = M Xn on C a (21.4) 

= 6 1 , on C u or M n n x + M*,n 2 =M 2n on C a (21.5) 

y/ x =y/ x on C u or M x n x +M X2 n 2 =M Xn on C a (21.6) 

y/ 2 = \j/~> on C u or Mfj n x +M 2 n 2 = M 2n on C a (21.7) 


where n x = cos (x x ,n) and n 2 = cos (x 2 ,n) are the components (direction cosines) of the unit outward 
normal vector to the cylindrical plate edges. 

The plate constitutive relations are derived by using Eqs. (2) and (3) with Eqs. (1 6)-(l 8), and 
integrating over the laminate thickness. The resulting constitutive relations of the new zigzag plate theory 
are expressed in matrix form as 

Nl r A B Olfe,; 

< M 6 11 D 0 < » (22) 

qJ 1_0 0 Gj[eJ 


15 



where expressions for the components of the stiffness matrices A = [/1 (/ ] ( 1,7 = 1 — 3), 
B = [#//] O’ = 1 — 3; y = 1 — 7) , D = [£E] (i,j = 1 - 7) , and G = [G jy ] (/,7 = 1-4) are given in 
Appendix A. 

Introducing Eqs. (22) into Eqs. (20) results in seven partial differential equilibrium equations in terms 
of seven kinematic variables, resulting in a 14 th -order theory. The equilibrium equations can be solved 
exactly or approximately depending on the complexity of the material lay-up, boundary conditions, and 
loading. In addition, because the highest partial derivative appearing in the strain measures in Eqs. (15) 
are first order, computationally efficient C°-continuous plate finite elements can be developed, thus 
enabling application of this zigzag theory in large-scale analyses of complex aerospace structures. 


6. Example problems and results 


To determine the accuracy of the present zigzag plate theory, analytic solutions for simply supported 
and cantilevered rectangular laminates are derived and detailed distributions of the displacements and 
stresses are examined. The rectangular laminates are defined on the domain 
x l e [0 ,a], x 2 e [0,6], z e [~h,h] . 

Example 1: A simply supported rectangular plate is subjected to the sinusoidal transverse pressure 
q(x l ,x 2 ) = q 0 sin(^x 1 /a)sin(;rx 2 lb ) . The simply supported boundary conditions are obtained from 
Eqs. (21). For cross-ply and uniaxial laminates, the kinematic and force boundary conditions along 
x, e [0,uf] are 


0 

II 

£ 

II 

II 

* 

II 

> 

(23.1) 

0 

II 

II 

II 

(23.2) 


and along x 2 e [0,6] , 


u = w = 6 X = y/\ = 0 


N 2 =M 2 = Ml = 0 


(23.3) 

(23.4) 


For this set of boundary conditions, the exact solutions are obtained by the following trigonometric 
expansions 


1 v = W sin (nx { / o)sin(^x 2 / 6) 


(24.1) 


u 


u' 

Ox 

> = < 

©i 





cos^Xj / a) sin(^x 2 / 6) 


(24.2) 
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> sm{nx x / a)cos{nx 2 / b ) 


(24.3) 


V 


' F ' 

o 2 

> = 5 

©2 

yi . 


.^2, 


where {[/, F, JF, 0 l5 0 2 , ¥ 2 } are the unknown magnitudes of the kinematic variables that are 

determined from satisfaction of the equilibrium equations. 

For antisymmetric angle -ply laminates, the kinematic and force boundary conditions along x 1 e [0,o] 
are 


u = w = 0 2 = y/ 2 = 0 

(25.1) 

o 

II 

ll 

sf 

ii 

CN 

(25.2) 

and along x 2 e [0,6] 


o 

II 

II 

II 

II 

> 

(25.3) 

N n =M 2 =M$ = 0 

(25.4) 


Thus, the trigonometric expansions that satisfy Eqs. (25) exactly differ from those in Eq. (24) only for the 
u and v variables, and they are given by 

u = U sin^Xj / a) cos (;rx 2 lb), v-V cos(srx 1 / a ) sin(;rx 2 / b) ^26) 

Example 2: A cantilevered rectangular plate clamped along a single edge (x, =0), free along the 
other edges, and subjected to the uniform transverse pressure q ( x, , x 2 ) = q o . The clamped boundary 
conditions along x t = 0 are 


u = v = w = 0 X = 0 2 = Y\ = W 2 = 0 (27 1) 

EJnlike the previous zigzag theories where the clamped boundary conditions result in erroneous solutions 
for transverse shear stresses and forces that vanish along the clamped edges, the solutions of the present 
theory do not possess such anomalies. For instance, along the clamped boundary at x x = 0 , the kinematic 
constraints in Eqs. (27.1), with the use of Eqs. (2.5), (13.2), and (18), give rise to the following 
transverse-shear stresses, r[ k z } (0, x 2 , z) , and force, Q x { 0, x 2 ) : 



(27.2) 
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where, in general, w l (0,x 2 ) ^ 0 . 


The traction-free boundary conditions along the edge x l = a are 
N x = N n =M X = M n = Q X = Mf = Mf 2 = 0 
and along the edge x 2 e [0,6] 


(28.1) 


N n =N 2 =M n =M 2 =Q 2 = Uf 2 =M* = 0 (28 .2) 

For both example problems, various laminates are considered with the emphasis on relatively thick 
laminated composite and sandwich plates having the span-to-thickness ratio a / 2h = b / 2h = 5 . The 
mechanical material properties are listed in Table 1, whereas in Table 2 the stacking sequences of the 
laminates are summarized. 


The example problems include: (1) a two-layer, cross-ply carbon-epoxy laminate, labeled as laminate 
A; (2) three variations of a three-layer sandwich laminate, laminates B, Bi, and B 2 , having uniaxial 
carbon-epoxy face sheets and a thick, closed-cell polyvinyl chloride (PVC) core. Flerein, PVC is 
represented as an isotropic material. Laminates B]_ and B 2 have progressively thinner face sheets; (3) a 
sandwich laminate with two- layered titanium face sheets and a thick titanium honeycomb core, laminate 
F; and (4) a five-layer, angle-ply sandwich laminate with carbon-epoxy face sheets and a thick PVC core, 
laminate G. 


For comparison purposes, additional analytic and finite element solutions were also obtained for the 
corresponding boundary-value problems using the three-dimensional elasticity theory [27, 28], FSDT, Di 
Sciuva theory [12], and MSC/MD-NASTRAN finite element code [29]. Note that application of FSDT 
generally requires the use of shear correction factors. For laminated composites, lamination-dependent 
shear correction factors have been shown to provide relatively accurate deflection predictions (e.g., refer 
to [30, 31]). Yet, such shear corrections fail to furnish substantial improvements for the normal strain and 

stress predictions. Presently, to establish a common framework reference for FSDT, k 2 =5/6, a shear 
correction factor that is appropriate for homogeneous plates, was used throughout. 

For the simply supported cross-ply and uniaxial laminates, the exact solutions for the FSDT and Di 
Sciuva zigzag theory were obtained using the trigonometric functions in Eqs. (25.3) (excluding those 
functions for y/ a ( a = 1,2) which do not appear in FSDT). For the angle-ply antisymmetric laminates, 

Di Sciuva theory permits only approximate solutions to be determined. Presently, the Rayleigh-Ritz 
method was employed, where the kinematic variables were approximated with suitable Gram- Schmidt 
polynomials [26] that satisfy the kinematic boundary conditions exactly, Eqs. (26.1) and (26.2) (refer to 
Appendix B for details on the Gram-Schmidt polynomials). Furthermore, for the simply supported cross- 
ply and angle-ply antisymmetric laminates, exact three-dimensional elasticity solutions were obtained 
using the solution procedures developed by Pagano [27] and Burton and Noor [28]. 

For the cantilevered laminates, approximate solutions corresponding to Di Sciuva and refined 
(present) zigzag theories were developed using the Rayleigh-Ritz method. Flere the displacement 
approximations are based on the Gram-Schmidt polynomials, using seven functions along the xi axis and 
five along the x 2 axis (see Appendix B). Furthermore, for the cantilevered plate (laminate B), a high- 
fidelity three-dimensional (3D) finite element solution was obtained by using MSC/MD-NASTRAN. The 
model is regularly discretized and is comprised of three elements through the thickness for each face 
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sheet, six elements through the core thickness, and sixty subdivisions along each span direction, with the 
total of 43,200 HEXA20, linear-strain elements. 

The numerical and graphical results that follow are labeled as: 

• 3D Elasticity (three-dimensional elasticity solutions using procedures developed by Pagano [27] 
for cross-ply laminates and by Noor and Burton [28] for angle-ply antisymmetric laminates). 

• FSDT (First-order Shear Deformation Theory, k 2 = 5 / 6 ). 

• Zigzag (D) (Di Sciuva [12]). 

• Zigzag (R) (present, refined zigzag theory). 

• 3D FEM (three-dimensional FEM solution using MSC/MD-NASTRAN [29]). 

• Zigzag (R-E) (transverse shear stresses obtained by way of integrating three-dimensional 
elasticity equilibrium equations, using the normal and in-plane shear stresses derived from the 
refined zigzag theory). 

Comparisons of the maximum deflection and maximum top-surface displacement are presented for the 
simply supported square laminates {a / 2h = b / 2h = 5) in Tables 3 and 4. These results demonstrate that 
both zigzag theories predict accurate plate displacements as compared to the three-dimensional elasticity 
solution, with the refined theory achieving slightly more accurate predictions. The laminate G case is 
somewhat pathological for Zigzag (D), because the solution in this case is only approximate and has not 
converged. The deflections predicted by FSDT are generally overly stiff; however, they are expected to 
improve substantially with the use of lamination-dependent shear correction factors. In Table 5, the range 
of applicability of the various theories is examined by comparing the maximum (central) plate deflection 
for a simply supported sandwich laminate B, where the solutions correspond to the span-to-thickness ratio 
in the range of 4 through 200. These results are also plotted in Figures 3(a) and 3(b). It is observed that 
both zigzag theories predict accurate deflections for all span-to-thickness ratios examined, whereas FSDT 
is overly stiff, especially in the thick regime, as expected. Even for a relatively thin laminate B 
(a /2h = 50), FSDT underestimates the maximum deflection by about 36%, which means that a much 
more significant shear correction factor would be required for FSDT for this type of laminate. 

For the simply supported laminates, normalized through-the-thickness distributions for the in-plane 
displacement, u l = (10 4 D 11 / q o a 4 ) u[ k) (0,a /2,z) , the normal stress, dj, = |(2/i)~ / /g 0 rr jx 

<Tjf \a/2,al2,z), and the transverse-shear stress, f lz = (2h/q 0 a)T^ ] (0,a /2,z) , are depicted in 
Figures 4-21. 

For laminate A - an asymmetric cross-ply carbon-epoxy composite - the u x and dj , quantities are 
accurately modeled by the FSDT, Zigzag (D), and Zigzag (R) theories, with Zigzag (D) producing 
slightly underestimated «] displacement near the top surface (Figures 4-6). The major result differences 

for this laminate correspond to the T lz distribution (Figure 6), where both FSDT and Zigzag (R) produce 
piecewise constant stresses, whereas Zigzag (D) gives a uniform (constant) distribution for r lz that is 

significantly less accurate than the predictions of the other two theories. In addition, Zigzag (R) theory 
provides a more accurate evaluation of the average transverse shear stress within each lamina than does 
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FSDT. Also, the r lz stress, obtained by equilibrium-based method. Zigzag (R-E), is best correlated with 
the three-dimensional elasticity solution (3D Elasticity). 

The w | , <J n , and Tj. results for sandwich panels B, Bj and B 2 that have constant-thickness carbon- 

epoxy face sheets and a PVC core are provided in Figures 7-15. The three laminates differ in terms of the 
face sheet thicknesses, starting with the thickest face sheets in laminate B and progressing to thinner face 
sheets in laminates B! and B 2 . This study evaluates the effect of face sheet thinness on the through-the- 
thickness displacement and stress variations. As evidenced from Figure 7, the u ] zigzag effect is very 
strong in laminate B, and is less pronounced as the face-sheet thickness diminishes, as shown Figures 10 
and 13 for laminates B! and B 2 , respectively. For laminate B, the top- and bottom-surface values of the u x 
displacement are significantly underestimated by FSDT that is only capable of a linear (average) 
distribution. Also, as shown in Figure 8, FSDT grossly underestimates the normal stress a u on the 
bounding surfaces where the greatest compression and tension occur. As the face sheets become thinner 
in relation to the core thickness, the FSDT predictions for u x and a n tend to improve, as evidenced from 

Figures 7, 8, 10, 11, 13, and 14. Fiowever, FSDT tends to overestimate the transverse shear stress, r lz , 

significantly, while underestimating this stress in the core (refer to Figures 9, 12, and 15). As seen from 
Figures 9, 12, and 15, Zigzag (D) provides accurate response predictions for all three sandwich panels, 
with the exception of the transverse shear stress; the T Xz distributions are accurate within the core but 

erroneous in the face sheets. By contrast, Zigzag (R) yields accurate solutions of all response quantities 
for the three sandwich panels examined. When the transverse shear stresses are evaluated from the 
constitutive relations, the theory provides the correct average values in the face sheets and in the core. 
Moreover, the equilibrium-based method. Zigzag (R-E), furnishes the superior transverse shear stresses 
that are virtually indistinguishable from those of three-dimensional elasticity. 

The sandwich study just presented illuminates some of the basic flaws in the classical sandwich- 
modeling assumptions, i.e., in-plane displacements are linear through-the-thickness and normal stresses 
are negligible in the core, and in-plane displacements are constant and transverse-shear stresses are 
negligible in the face sheets (e.g., refer to [32-34]). The hypothesis concerning the normal stresses is 
generally correct in the range of face-sheet thicknesses examined, as seen in Figures 8, 11, and 14; 
whereas the assumptions concerning the in-plane displacements and the transverse shear stresses are only 
confirmed for the sandwich panels with very thin face sheets. As seen in Figure 7, which depicts the 
laminate B results, the in-plane displacement exhibits a high gradient through the face sheets, thus 
departing appreciably from the classical, membrane assumption of uniform displacement. Furthermore, 
the average transverse shear stress in the face sheets is higher than in the core (see Figure 9), again 
contrasting with the classical assumption of negligible transverse shear stress in the face sheets. 

The displacement and stress results for laminate F - a cross-ply sandwich with a lower degree of 
anisotropy than laminate B - are depicted in Figures 16-18. For this laminate, the zigzag effect of the in- 
plane displacement is somewhat less pronounced than for laminate B. As in the previous example, Zigzag 
(R) yields highly accurate predictions of all response quantities. It is seen that FSDT provides over 
estimated values for the in-plane displacement on the bounding surfaces. Also, it is evident from Figure 
18(a) that the integral of the shear stress over the thickness computed using the Zigzag (D) transverse 
shear stress would result in a significantly greater value than the exact shear force. 

Figures 19-21 demonstrate the results for laminate G - an angle-ply antisymmetric sandwich plate 
with multilayered face sheets. This is a highly challenging test case for any lamination theory. For this 
lamination, only an approximate solution can be obtained for Zigzag (D), requiring a large number of 
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suitable shape functions to achieve a converged solution. Using the Rayleigh-Ritz method with the Gram- 
Schmidt polynomials approximating the kinematic variables (see the details of the approximating 
functions in Appendix B), a relatively inaccurate solution was obtained. Consequently, the Zigzag (D) 
results for the in-plane displacement (Figure 19) and normal stress (Figure 20) are somewhat erroneous. 
On the other hand. Zigzag (R) enables an exact solution to be obtained for this problem, once again 
yielding highly accurate predictions of the response quantities. As in the previous examples, FSDT 
models the in-plane displacement response only in an average sense (a linear distribution through the 
thickness), leading to a significant underestimation (a non conservative prediction) of the normal stress in 
the face sheets. 

To examine the effect of clamped boundary conditions, a square cantilevered laminate B under a 
uniform transverse pressure was examined. Table 6 summarizes the maximum deflection calculated using 
the three different theories and a three-dimensional finite element solution that serves as a reference. The 
FSDT deflection is underestimated by an order of magnitude. The two zigzag theories give accurate 
results, with Zigzag (R) producing a somewhat superior deflection prediction. 

For the cantilevered sandwich laminate B, normalized through-the-thickness distributions of the in- 
plane displacement, wj = (10 4 D n / q Q a 4 ) u[ k \a,a /2,z) , the normal stress, <t u =|(2 h)' /q 0 a 2 ^jx 

CTj f ’ (a / 5, a / 2, z) , and the transverse-shear stress, rj 7 = (2 h/q 0 a)x r[ k \a / 5,a /2,z) , are provided in 

Figures 22-24. The stresses were computed near the clamped edge (x, =a/5,x 2 = a/ 2) to allow for 

proper comparisons with the accurate stresses obtained from the 3D FEM analysis. For this problem, both 
zigzag theories produce accurate results; however, application of Zigzag (R) resulted in superior 
predictions of transverse shear stresses. Finally, the normalized transverse shear force, 

Q x = — [ t\\ 1 tlz , evaluated at x 2 = a 12 , is plotted versus the normalized axial coordinate, x, la , as 
Vo a ' h 

shown in Figure 25. For this problem, both FSDT and Zigzag (R) predict the correct linear distribution, 
yielding a maximum value at the clamped edge and vanishing at the free edge. This contrasts with an 
erroneous Zigzag (D) solution that varies in a non-linear manner across the span and which vanishes at 
the clamped edge. 


7. Conclusions 


A refined zigzag theory has been developed for laminated-composite and sandwich plates that exhibit 
a high degree of transverse shear flexibility, anisotropy, and heterogeneity. In this refined theory, a first- 
order shear-deformation theory is used as a baseline for the kinematic assumptions with a set of novel 
piecewise-continuous zigzag displacements added to the in-plane displacement components. The resulting 
kinematic field is independent of the number of material layers, and the zigzag displacements are defined 
by requiring only partial lamina- interface continuity requirements of transverse shear stresses. The force 
equilibrium equations, boundary conditions, and strain-displacement relations are completely consistent 
with respect to the virtual work principle, and transverse-shear correction factors are not required. The 
refined zigzag theory is better suited for engineering practice than previous similar theories because of its 
relative simplicity and its ability to model accurately the transverse shear and in-plane deformations of the 
individual laminae in a physically realistic manner. Unlike other similar theories, meaningful in-plane and 
transverse shear stresses are obtained directly from the constitutive equations, in a theoretically consistent 
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manner. The new theory is devoid of a major shortcoming of other similar theories; that is, the new theory 
enables accurate modeling of clamped boundary conditions. 

Results for several example problems have been presented that highlight the superior predictive 
capability attainable with the present theory and its ability to model correctly clamped boundary 
conditions. The critical quantitative assessment of the new theory, that included analyses of highly 
heterogeneous sandwich laminates in bending, revealed that this refined zigzag theory is more accurate 
than previous similar theories. 

An additional and important benefit is that the new zigzag theory lends itself well for finite element 
approximations. In particular, the theory is perfectly suited for the development of computationally 
efficient, C°-continuous finite elements. Because of a wide applicability range that includes laminated- 
composite and sandwich structures, such finite elements would be highly desirable for large-scale 
analyses and design studies of high-performance aerospace vehicles. 
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Appendix A. Plate stiffness coefficients 


The stiffness coefficients in Eq. (24) are computed from the following expressions using the 
constitutive coefficients given in Eq. (3). 

(1) Matrix A = [ A if - ] (i,j = 1-3) , symmetric, (3x3): 

A=[ Cdz (A.l) 

J-h 

(2) Matrix B = [ B t] ] (/ = 1 - 3; / = 1 - 7) , non-symmetric, (3x7) : 

B = I* CB^ dz (A.2) 


(3) Matrix D = [ZG] (i,j = 1 - 7) , symmetric, (7x7) : 

ds 1>; C b^ 

(4) Matrix G = [G’ ;/ ] (/, j = 1 - 4) , symmetric, (4x4) : 

G = Bp Q Bp Jz 

where 
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Appendix B. Gram-Schmidt polynomials 


The Gram-Schmidt orthogonal polynomials [26] are widely used as shape functions for application to 
the Rayleigh-Ritz method. These functions can accommodate various kinematic boundary conditions and 
yield a diagonal mass matrix. 

The procedure for constructing such polynomials is initially established for one-dimensional domains, 
whereas for two-dimensional applications, simple products of one-dimensional functions are used. 

The general expression of a one-dimensional Gram-Schmidt polynomial %(s), where .ye [0, 1 ] is a 
dimensionless axial coordinate, is 

X p+ i( s )={ s -A p )z p (s)-B p Z p -i(s) (B. 1 ) 

where the coefficients A p and B p are computed from the expressions 

\s Zp \s)ds \ y 2 {s)ds (B.2) 

A = ll R =_l2 

p ~ r 1 7 ’ d p _ r i , 

j 0 Z P (s)ds J o Z P -i (s)ds 

The process of building these functions initiates with the first two polynomials 


Xo( s ) = 0 

Zi(s) = M s f 2 


(B.3) 


where 


A (s) = s 
jU 2 (s) = \-s 


(B.4) 


The three exponent values O, = (0,1,2) correspond, respectively, to the end condition of the function: 0 

if the function does not vanish, 1 if the function vanishes, however, its first derivative does not vanish, 
and 2 if the function and its first derivative vanish. 

The orthogonality of the Gram-Schmidt polynomials is manifested by the relationship 

X P (s)z 9 (s)ds = c pq 8 pq (B-5) 
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where c pq is a coefficient and 8 pq is the Kronecker’s delta. 

For the cantilevered plate examined in Section 6 , the Gram-Schmidt polynomials employed for the 
kinematic variables have the form 

M 1 V 

/ (*1 > *2 ) = Z Z fijZl ^ )Zj ( x 2 ) ( B -6) 

i = 1 j = 1 

where M=1 , N= 5, and are the u nk nown amplitude coefficients to be determined from the Rayleigh- 

Ritz analysis. 


Zi(x l ) = x l 


Q, =1 
Q 2 =0 


I 2 ( x i ) = x l-^ x i 

, , 3 4 2 2 

Z,(x } ) = x L --x L + -x 1 

M 4 15 3 ^ 15 2 5 

y d (x,) = x, x, + — x, x, 

— 8 — 14 — 28- 

/ \ 5 12 4 _ 3 2 2 1 
2 T 5 (x L ) = x L -yX L + 2 x L --x L +— ^ 
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(B.7) 


Q, =0 

Q 2 =0 


Z t ( x l ) = l 
Z 2 (x 2 ) = x 2 
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j 3 (x 2 ) = x 2 -x 2 +^ 
— — — 6 


/ ^ 3 3 2 3 

J 4 (x 2 ) = x 2 --x 1 + -x 2 - — 
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Z 5 (x 2 ) = x 2 - 2 x 2 + -x 2 2 - ^x 2 + 


1 _ 
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— X 2 H 

— — — 7 — 7 — 70 
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(B. 8 ) 
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(B.9) 


x, = xj a 
x 2 = x 2 / b 

For the simply supported angle-ply antisymmetric laminates considered in Section 6, the Di Sciuva 
zigzag theory solutions utilizing the Rayleigh-Ritz method employ the Gram-Schmidt polynomials 
corresponding to M=N=5. Because the kinematic variables have different boundary conditions, Eq. (26.1- 
a) and Eq. (26.2-a), the Gram-Schmidt polynomials are chosen according to the table below. 


Table B. Parameters for Gram-Schmidt polynomials. 


Kinematic variable 

zOi) 

Zi(x 2 ) 


1^2 


f} 2 

u 

1 

1 

0 

0 

V 

0 

0 

1 

1 

w 

1 

1 

1 

1 

¥1 

0 

0 

1 

1 

¥2 

1 

1 

0 

0 


For example, the Xi-direction functions for the transverse displacement, w, have the form, 


Q, = 1 

Q 2 = l 


Zi(x 1 ) = xl~x 1 

, , 3 3 2.1 

Z 2 (x i ) = x L --x L + -x i 
f \ 4 o 3 17 2 3 

Z,(^ 1 ) = x l - 2 x l + —x l -—x 1 

/x 5 5 4 13 3 3 2 1 

j 4 (^) = Xl - -Xl + — jq - -Xl + — Xf 

, , 6 , 5 37 4 19 3 13 2 1 

Zs(x i) = ^L-3^ L + YYf L _ YT— 33 — - 33 - 


(B.10) 


Similar functions are constructed for the x 2 -direction. 
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Table 1. Mechanical material properties. 


Orthotropic materials 

Isotropic materials 

Lamina 

material 

C 

H 

Lamina 

material 

P 

T 

Carbon-Epoxy 

unidirectional 

composite 

Titanium 

honeycomb 

core 

PVC core 

Titanium 

E[ k) [GPa] 

1.579 x 10 2 

1.915 x 10' 1 

E (k) 

[GPa] 

1.040 x 10' 1 

1.041 x 10 2 

E (k) [GPa] 

9.584 

1.915 x 10' 1 

E™ [GPa] 

9.584 

1.915 

v (k) 

v \2 

0.32 

0.658 xlO' 2 


0.3 

0.31 

v (k) 

M3 

0.32 

0.643 xlO' 6 

y(t) 

k 23 

0.49 

0.643 xlO' 6 

G (k) [GPa] 

5.930 

4.227 xlO' 8 


G, (k) [GPa] 

5.930 

5.651 x 10' 1 

G™ [GPa] 

3.227 

1.248 
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Table 2. Laminate stacking sequences (lamina sequence is in the positive z direction). 


Laminate 

Normalized lamina 
thickness, h' k> /h 

Lamina 

materials 

Lamina orientation 
[°] 

A 

Cross-ply 

composite 

(0.5/0. 5) 

(C/C) 

(0/90) 

B 

Uniaxial 

sandwich 

(0.1/0. 8/0.1) 

(C/P/C) 

(0/0/0) 

Bi 

Uniaxial 

sandwich 

(0.025/0.95/0.025) 

(C/P/C) 

(0/0/0) 

b 2 

Uniaxial 

sandwich 

(0.0025/0.995/0.0025) 

(C/P/C) 

(0/0/0) 

F 

Uniaxial 

sandwich 

(0.1/0. 8/0.1) 

(T/H/T) 

(0/0/0) 

G 

Angle-ply 

sandwich 

(0.05/0.05/0.8/0.05/0.05) 

(C/C/P/C/C) 

(30/-45/0/45/-30) 


Table 3. Normalized maximum (central) deflection, w = (10 2 D U / q 0 a 4 )w(a /2,b/2), 
for simply supported laminates subjected to sinusoidal transverse pressure loading. 


Laminate 

Normalization 

factor, 

10 2 D n /q Q a 4 

3D 

Elasticity 

FSDT 

Zigzag (D) 

Zigzag (R) 

A 

8.147 x 10' 2 

1.228 

1.278 

1.170 

1.219 

B 

7.502 xlO' 2 

29.761 

2.731 

29.769 

29.785 

Bi 

2.201 x 10' 2 

11.645 

2.819 

11.693 

11.694 

b 2 

2.402 xlO' 3 

2.080 

1.728 

2.103 

2.103 

F 

5.444x1 O' 2 

1.331 

0.389 

1.332 

1.333 

G 

3.55 lxlO' 2 

14.124 

1.055 

12.734 

14.105 


30 





Table 4. Normalized maximum (top surface) displacement along x 1 axis u l = (10 3 Z) n / q Q a 4 ) x 
u\ N \a,b 1 2,h) of simply supported laminates subjected to sinusoidal transverse pressure loading. 


Laminate 

Normalization 

factor, 

10 3 D n /q Q a 4 

3D 

Elasticity 

FSDT 

Zigzag (D) 

Zigzag (R) 

A 

8.147 x 10' 1 

4.233 

4.152 

3.855 

4.251 

B 

7.502 xlO' 1 

9.977 

2.156 

9.945 

9.897 

B1 

2.201 x 10' 1 

2.008 

2.099 

2.030 

2.032 

B2 

2.402 xlO' 2 

1.809 

1.993 

1.851 

1.851 

F 

5.444x1 O' 1 

0.643 

0.796 

0.646 

0.649 

G 

3.551 x 10' 1 

3.908 

0.704 

1.295 

3.845 


Table 5. Normalized maximum (central) deflection, w = (10 2 Z) 11 / q 0 a 4 )w(a / 2,b / 2) , for simply 

supported laminate B subjected to sinusoidal transverse pressure loading and corresponding to various 

span-to-thickness ratios. 


Span-to- 

thickness 

ratio 

Normalization 

factor, 

10 2 D n /q 0 a 4 

3D 

Elasticity 

FSDT 

Zigzag (D) 

Zigzag (R) 

4 

1.832 xlO' 1 

42.420 

3.739 

42.124 

42.189 

10 

4.668 xlO' 3 

9.734 

1.321 

9.738 

9.739 

20 

2.93 lxlO’ 4 

3.487 

0.948 

3.489 

3.490 

50 

7.502 xlO' 6 

1.305 

0.841 

1.305 

1.305 

100 

4.688 xlO' 7 

0.945 

0.826 

0.945 

0.945 

200 

2.93 lxlO' 8 

0.852 

0.822 

0.852 

0.852 
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Table 6. Normalized maximum (free-edge) deflection, w = (10 2 D n / q 0 a 4 )w(a,b/2 ) , for cantilevered 
laminate B subjected to uniform transverse pressure loading. 


Normalization 

factor, 

10 2 D n /q 0 a 4 

3D FEM 

FSDT 

Zigzag (D) 

Zigzag (R) 

7.502 xlO' 2 

246.778 

25.351 

244.077 

245.615 
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Span-to-thickness ratio, a/2h 

Figure 3 (a). Normalized central deflection vs. span-to-thickness ratio for simply supported laminate B 

subjected to sinusoidal transverse pressure. 



Span-to-thickness ratio, a/2h 

Figure 3 (b). Normalized central deflection vs. span-to-thickness ratio for simply supported laminate B 
subjected to sinusoidal transverse pressure (zoomed view). 
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Displacement, u x {0,cil 2,z) 

Figure 4. Normalized in-plane displacement for simply supported laminate A subjected to sinusoidal 

transverse pressure. 



Figure 5. Normalized normal stress for simply supported laminate A subjected to sinusoidal transverse 

pressure. 
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Figure 6. Normalized transverse shear stress for simply supported laminate A subjected to sinusoidal 

transverse pressure. 



Figure 7. Normalized in-plane displacement for simply supported laminate B subjected to sinusoidal 

transverse pressure. 
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Stress, cr n (a/2,a/2,z) 

Figure 8. Normalized normal stress for simply supported laminate B subjected to sinusoidal transverse 

pressure. 



Stress, r lz (0,a/2,z) 

Figure 9. Normalized transverse shear stress for simply supported laminate B subjected to sinusoidal 

transverse pressure. 
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Displacement, z^(0,a/2,z) 

Figure 10. Normalized in-plane displacement for simply supported laminate Bi subjected to sinusoidal 

transverse pressure. 



Figure 11 (a). Normalized normal stress for simply supported laminate Bi subjected to sinusoidal 

transverse pressure. 
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Stress, cr n (a/2,a/2,z) 

Figure 1 1 (b). Normalized normal stress near the top layer for simply supported laminate B , subjected to 

sinusoidal transverse pressure. 



Figure 12 (a). Normalized transverse shear stress for simply supported laminate B , subjected to 

sinusoidal transverse pressure. 
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Stress, z"| z (0,a/2,z) 

Figure 12 (b). Normalized transverse shear stress near the top layer for simply supported laminate B 

subjected to sinusoidal transverse pressure. 



Displacement, W[(0,a/2,z) 

Figure 13. Normalized in-plane displacement for simply supported laminate B2 subjected to sinusoidal 

transverse pressure. 


39 


0.5 


o3 

O 

O 

o 

C/3 

C/3 

<D 

2 

H 


-0.5 


-1 


(C, 0 U ) 


r 


-o- 


(P) 


— 3D Elasticity 
--FSDT 
— Zigzag (D) 
■^Zigzag (R) 


(C, o°) 


J 


-30 


-20 


-10 


10 20 


30 


Stress, cr n (a/2,a/2,z ) 

Figure 14 (a). Normalized normal stress for simply supported laminate B 2 subjected to sinusoidal 

transverse pressure. 



Stress, <r n (a/2,a/2,z) 

Figure 14 (b). Normalized normal stress near the top layer for simply supported laminate B 2 subjected 

to sinusoidal transverse pressure. 
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Figure 15 (a). Normalized transverse shear stress for simply supported laminate B 2 subjected to sinusoidal 

transverse pressure. 



Stress, r lz (0,a/2,z) 

Figure 15 (b). Normalized transverse shear stress near the top layer for simply supported laminate B 2 

subjected to sinusoidal transverse pressure. 
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Displacement, « 1 (0,a/2,z) 

Figure 16. Normalized in-plane displacement for simply supported laminate F subjected to sinusoidal 

transverse pressure. 



Stress, cr u (a/2,a/2,z) 

Figure 17. Normalized normal stress for simply supported laminate F subjected to sinusoidal transverse 

pressure. 
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Figure 18 (a). Normalized transverse shear stress for simply supported laminate F subjected to sinusoidal 

transverse pressure. 



Stress, zj z (0,a/2,z) 

Figure 18 (b). Normalized transverse shear stress near the top layer for simply supported laminate F 

subjected to sinusoidal transverse pressure. 
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Displacement, W[(0,a/2,z) 


Figure 19. Normalized in-plane displacement for simply supported laminate G subjected to sinusoidal 

transverse pressure. 
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Stress, <r n (a / 2,a / 2 9 z) 

Figure 20. Normalized normal stress for simply supported laminate G subjected to sinusoidal transverse 

pressure. 
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-0.5 0 0.5 1 1.5 

Stress, rj,(0,a/2,z) 

Figure 21 (a). Normalized transverse shear stress for simply supported laminate G subjected to 

sinusoidal transverse pressure. 



Stress, rj z (0,a/2,z) 

Figure 21 (b). Normalized transverse shear stress near the top layer for simply supported laminate 
G subjected to sinusoidal transverse pressure. 
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Figure 22. Normalized in-plane displacement for cantilevered laminate B subjected to uniform transverse 

pressure. 



Figure 23. Normalized normal stress for cantilevered laminate B subjected to uniform transverse pressure. 
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Stress, T l: (a/ 5,a/2,z) 

Figure 24 (a). Normalized transverse shear stress for cantilevered laminate B subjected to uniform 

transverse pressure. 



Stress, r lr (a/5,a/2,z) 

Figure 24 (b). Normalized transverse shear stress for cantilevered laminate B subjected to uniform 

transverse pressure. 
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Axial coordinate, x x t a 

Figure 25. Normalized transverse shear force along centerline of cantilevered laminate B subjected to 

uniform transverse pressure. 
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